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Abstract 

The recently observed L'*j(2317)^ meson is explained as a scalar cs system which ap- 
pears as a bound state pole, below threshold, in the DK scattering amplitude. The standard 
cs charmed scalar DgQ is found at about 2.9 GeV, with a width of some 150 MeV. 

For the scattering length of DK in S wave we predict 5 it 1 GeV~^. 



1 Introduction 

The 23 17)"*" state has, since its discovery [1] by the BABAR collaboration (see also Refs. [2,3]), 
raised several interesting questions on its nature [4-15]. However, the light-scalar- meson hypoth- 
esis from unitarization [16] has made important progress with the discovery of the D*j(2317)"'" 



1 



state (most likely J — 0). Not often it happens that in experiment a state is discovered which 
seems difficult to handle within other theories [17-21], but can be explained by the unitarization 
hypothesis without extra ingredients. 

Of course, this happened before with the light scalar mesons, which can easily be handled 
within a unitarized model for all non-exotic mesons [22] , without making any distinction between 
the heavier scalar mesons and the light ones. But seventeen years after the publication of these 
results [16], their implications are still not fully understood, nor applied in data analysis. 

In the specific case of the D*j{2317)~^ state, just by applying the model parameters established 
twenty years ago [23] , the unitarization scheme results in a bound state below the DK threshold, 
for a cs confinement spectrum with ground state at 2.545 GeV, almost 200 MeV above the DK 
threshold. This achievement is entirely similar to the extra states obtained by the same model 
for the light scalars. No additional interactions are necessary to explain the scalar mesons. 

2 Unitarization 

In the unitarization scheme, we describe meson-meson scattering by the wave equation 

(f ) = {E- Hf)-' Vt {E - H,)-' Vt ^Pf (f ) . (1) 

The operators Hf and respectively describe the dynamics of two free mesons and of a confined 
quark-antiquark system. The latter is supposed to generate an infinity of confinement states. 
Transitions from confined quarks to free mesons, mediated by quark-antiquark-pair creation, is 
described by Vt- 

By comparison of Eq. (1) with the usual expressions for scattering wave equations, we conclude 
that the generalized potential V for meson-meson scattering is given by 

V ^ Vt{E- H,)-' Vt . (2) 

With the definitions given in Appendix A, in particular formulae (7), (17) and Table 1, we 
obtain for the partial- wave scattering phase shift generated by wave equation (1) the result 

n=0 ^yP) ~ ^ne. 

The expression (3), which is an exact solution of Eq. (1), is very compact and elegant. But the 
summations over the radial excitations of the confined system do not converge fast. Hence, for 
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the numerical calculation of scattering quantities it is not suitable. A more practical method 
has been described in Ref. [24]. Nevertheless, because of its close contact with the physical 
situation, expression (3) is extremely comprehensible, and useful as a starting point for making 
approximations . 

For example, when we are studying a limited range of energies, we may restrict the summation 
to only the few terms near in energy, characterized by Eni^ in Eq. (3), and approximate the 
remaining terms by a simple function of p, the relative meson-meson linear momentum. Notice, 
however, that the resulting expression is only of the type of a Breit-Wigner (BW) expansion when 
the strength of the transition, characterized by A, is small. For larger values of A, expression (3) 
behaves very differently from BW expansions. 

Another, in practice very useful approximation is inspired by the radial dependence of the 
transition potential Vj, which is depicted in Ref. [25] for meson-meson scattering in the cases 

— 0~, — 1~, and — O"*", We find there that Vt is peaked at short distances. If we, 
moreover, take into account that eigenfunctions of the confinement Hamiltonian He must also 
be of short range, then we may just define a transition radius a and approximate the spatial 
integrals of formula (17) by choosing a spherical delta shell U{r) = S{r — a). In this case we end 
up with the expression 



2a^X^^ip jiipa) n(,{pa) ^ 



C0tg((5^(p)) oo i-r- ^ ^ 2 ■ (4) 

n=0 ^yP) ~ ^nic 

This way we have pulled the p dependence outside the infinite summation over the radial con- 
finement spectrum, which makes it much easier to handle truncations, as now the rest term docs 
not depend on p and may thus be chosen constant. Our recent detailed analyses of scalar mesons 
have actually been performed in the approximation of formula (4) . 
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Figure 1: Theoretical cross section cr as a function of the strength of the transition potential, 
characterized by A, for Kti elastic S'-wave scattering, coupled to a confinement spectrum with 
ground state at 1.389 GeV and level spacing of 380 MeV (model parameters taken from Ref. [23]). 

In Fig. (1) we show how the model's cross section varies with A in the approximation given 
in Eq. (4), using harmonic-oscillator confinement. The first two plots (A = 0.2 and 0.3) show 
separate peaks, not very different from a sum of BW resonances. For larger values of A, the cross 
section gets more and more distorted, certainly not suited for a description in terms of separate 
or overlapping BW resonances. The physical value for the strength of the transition potential 
comes in our model at around A = 0.75, for which value we find a broad structure peaking at 
some 830 MeV, representing the recently confirmed [16,26] Kq(800) (or k) resonance. Note that 
the latter structure is already showing up for smaller coupling, and is clearly not related to the 
ground state of the = 0^ sd confinement spectrum, at 1.389 GeV for the model parameters 
from Ref. [23]. 

The plots for A = 0.75 and A = 1.0 in Fig. 1 show structures in the scattering cross section that 
are certainly not well approximated by sums over BW resonances. Nevertheless, expression (3) 
has several features in common with BW resonances. 

A BW resonance is characterized by a singularity in the lower-half of the complex-energy 
plane (second Riemann sheet for single-channel scattering) . The real part Eq and the imaginary 
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part r/2 of the pole are about equal to the central resonance position and half the width of the 
resonance, respectively. 

In our model, as we have the disposal of an analytic expression, given in Eq. (3), we can 
determine the pole positions exactly. For small values of A, the relation between the central peak 
and width of a structure in the cross section can rather well be read from the real and imaginary 
part of the pole. However, for larger values of A, poles of the scattering amplitude are not always 
close to the corresponding peaking structure in the cross section. 

Nevertheless, poles are very useful, because they constitute a clean tool to relate phenomena 
of different hadronic systems [16,27-30]. For example, only through the hypothetical movements 
of the poles in the complex-energy plane we can show that the /o(600), /o(980), ^0(800) [26] 
and ao(980) are manifestations of the same phenomenon, but for different flavors [31]. 

3 The scalar mesons 

The scalar mesons represent the most noteworthy result [16] of the unitarized model given in 
formula (1). In the recent past, scalar mesons have been the subject of many investigations, 
in a variety of approaches. Strong and electroweak interaction properties of scalar mesons were 
studied in Refs. [32-34], quark and glue contents of the isoscalars /o in Refs. [35-50], the ft;(800) 
in Refs. [31,32,51-53], the isovectors Oq in Ref. [54,55], dynamically generated scalar mesons 
in Ref. [56], analogy with the Higgs boson in Refs. [57-59], effects of scalar mesons in decay 
processes and coupling to meson-meson scattering in Refs. [60-62], ao//o mixing in Ref. [63-66], 
mixing of hght and heavy scalar mesons in Refs. [67,68], and properties of scalar mesons at finite 
temperatures in Ref. [69]. 

Our understanding of the scalar mesons is based on studying the corresponding pole positions, 
and their movements as a function of the coupling A. A pole which stems from a confinement 
state moves towards the energy eigenvalue of that state when A is changed to smaller values. 
However, we observe that besides the poles stemming directly from the confinement spectrum, 
there exist poles in the scattering amplitude which disappear into the background for vanishing 
A. These somehow originate in the bulk of the sum over the radial excitations in expression (3), 
and are not related to any state in particular. 

Moreover, we observe that the lightest i = 1 qq scalar mesons can only reside somewhere 
between 1.3 and 1.5 GeV in the confinement spectrum, with the model parameters of Ref. [23]. 
The latter have been adjusted to bottomonium, charmonium, the light vector states, and to 
P-wave meson-meson scattering data, whenever available in the early eighties. 

Finally, we observe that the extra poles in the scattering amplitude come close enough to the 
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real-energy axis only in the case of S-weive scattering. There, they show up somewhere close to 
the lowest threshold, for energies below the ground states of the confinement spectrum. 

For the physical value of the model parameter A, a specific value has been selected by Nature. 
But in our model we can vary A. This way, we can find out where poles are moving for smaller 
or larger values of the coupling. 



4 The (2317)+ state 



sJ 

For DK elastic S-wave scattering, we adopt the same expression that has been used for Kn as 
described in Ref. [31]. We only have to substitute the confinement masses of ds by those for cs, 
thereby using the parameters of Ref. [23]. For the ground state of the = 0^ cs confinement 
spectrum, the mass comes out at 2.545 GeV. The first radial excitation is then 380 MeV (the 
model's level spacing) higher in mass. The summation over the higher radial excitations in Eq. (4) 
we approximate by a constant, normalized to 1. The relative couplings of the ground state and 
the first radial excitation to the DK channel are 1.0 and 0.2, respectively [31]. Thus we obtain 

|2 



E 



|^n,l(«) 



1.0 



+ 



0.2 



- 1 GeV- 



(5) 



E{p) - K,i E{p) - 2.545 E{p) - 2.925 

For expression (4) with the substitution (5), and furthermore for the physical value of A (meaning 
the value fitting Kn elastic 5"- wave scattering from threshold up to 1.6 GeV [31]), we depict the 
resulting DK elastic 5"- wave cross section in Fig. 2. The corresponding pole positions in the 
complex-energy plane for the D*q ground state and first radial excitation are shown in Ref. [6]. 
These are found at 2.28 GeV and (2.78 - i 0.093) GeV, respectively . 
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Figure 2: Theoretical cross section a for DK S-wave elastic scattering (model parameters taken 
from [23]). 



Notice from Fig. 2 that in the region around 2.55 GeV, where quark models predicted — 0"*" 
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cs states [17,21], we find no structure in the cross section, nor poles in the scattering amphtude. 
On the contrary, employing confinement plus unitarization, we find structures in the cross section 
some 250 MeV below and 400 MeV above that value, which explains the observation of the 
(probable) D*o ground state below the DK threshold, and predicts for the first radial excitation 
a resonance somewhere between 2.8 and 3.0 GeV, with a width of some 100-200 MeV. Based on 
Fig. 2, wc may also predict the scattering length for DK elastic S'-wave scattering. Keeping in 
mind that the = O"*" cs ground state comes out a little bit too low in mass with our value of 
A, we estimate the theoretical error by varying A. For a conservative choice, we obtain 

scattering length DK in S wave = 5 ± 1 GeV"^ . (6) 

5 Conclusions 

The two most important aspects of strong interactions for hadrons are confinement and OZI- 
allowed real or virtual decay, which effectively describe gluonic interactions and quark-pair cre- 
ation/annihilation, respectively. These are the only ingredients entering the model defined in 
Eq. (1), which, with a fixed small set of parameters, can reasonably describe hadronic phenom- 
ena ranging in energy from the light scalars up to bottomonium. 

Here, we have discussed the model's results for the — 0"*" cs states coupled to the DK 
channel. For the DK 5*- wave scattering length we predict the value 5 ± 1 GeV~^. 

A The T matrix for meson-meson scattering 

First, we introduce the operators He, describing the confinement dynamics in the interaction 
region, Hf, representing the dynamics of the scattered particles at large distances, and Vt, which 
stands for the transitions between these two sectors. For an arbitrary spherically symmetric 
confinement potential Vc, we define these operators in configuration space by 

and Vt = XU{r) . 

The various mass parameters of Eq. (7) are defined in Table 1. The transition potential Vt, 
which provides the communication between the confined channel and the scattering channel, is 
in agreement with the potentials [25] that may describe the breaking of the color string. 
Before performing the calculations to all orders, we write the transition operator as 

T = (I -VGfr'V (8) 
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symbol 


definition 


777,„ (fn^) 
It tq yiiLqj 

1^ 


ronstitiipnt ( anti- iniiark mass 

U X U LAV-' XX U \ CAiXX U X J VJ LACAiX XV XXXCAil_J>-/ 

reduced mass in confinement channel 
meson masses 

reduced mass in scattering channel 



Table 1: Mass parameters used in Eq. (7). 



= { 1 + VGf + VGfVGf + VGfVGfVGf + ■■■}¥= ^ {VGfy V , 

N=0 

where Gf represents the free boson propagator associated with the Hamiltonian Hf in Eq. (7), 
given by 



Gf Ik ,k;z 



k 



k 



(9) 



The matrix elements of the T operator follow from Eq. (8) by 



/oo 
d'k y: {p\ (VGf) 

N=0 



N 



k) Ik 



V 



p 



(10) 



p 



(11) 



For the generalized potential (2), the Born term takes the form 

V{p,p') = {p\v\p') = {p\Vt {E{p)-H,)-' Vt 
The total center-of-mass energy E and the linear momentum p are, through Eq. (7), related by 

(12) 



->2 
P 



E{p) = ^ + Ml + M2 . 

We denote the properly normalized eigensolutions of the operator He in Eq. (7), corresponding 
to the energy eigenvalue E^t^, by 

n = 0, 1, 2, 

(f |n4mc) = ytcmcif) J^niXr) , with 



. , ^, ... 
0. 1. 2. . . . 



(13) 



m = -4, .... +4 

So, by letting the self-adjoint operator He act to the left in Eq. (11), we find 

(P I y \p) = E (^^1 Vt \ntjne) {neerriel {E{p) - H,)-' Vt 

n£cmc 



P 



2^ {p\Vt , — - — Vt 



p 



The individual matrix elements in the sum of formula (14) result in 

(n^cmd Vt \k) = A {i)' Jneik) y;^ (k) , 
8 



(14) 



(15) 



where, in order to simplify the expressions, we define the quantities 

^ J r^dr ji ipr) U (r) J^*^^ (r) 



Jniip) 



r'^dr Ui ipr) U (r) JT* (r) 



(16) 



At this stage we should pay some attention to the quantum numbers of the two coupled 
systems, since the operator Vt mediates between systems with different quantum numbers. Prom 
parity conservation we deduce that the relative orbital angular momenta £ of a meson-meson pair 
and £c of the confined quark-antiquark pair differ by at least one unit. Also the total intrinsic 
spins are different for the two systems, when C parity is taken into account. The conservation of 
J, P, C, isospin, fiavor, and color can be formulated in a consistent way [25,70,71], but we do 
not intend to go into details here. We just should recall that there exists a relation between (£c, 
rric) and (£, m), and moreover that the overlap integral expressing this relation is absorbed into 
A. 

For the full Born term (14), we then obtain 



{p\V 



P' 



(17) 



Similarly, one may determine the A^+l-st-order term (p | [VGf)^^^ k) of Eq. (10), yielding 



N+l 



P 



An 



E (2^ + 1) 



TV 



X 



X Pe{p-p') Yl 



E{p)-En't, 2flfZ-p' 



/2 



(18) 



From Eqs. (17) and (18), we find for the matrix elements of the product (VGf)^ V the expression 



d^k {p\ {VGf) 



N 



k) Ik 



V 



p ) = 



nic 



X Pe{p-p') J2 



N 



X 



(19) 



E{p) - E^,e. 

The terms labelled with N in Eq. (19) can be summed up, resulting in a compact final formula 
for the matrix elements of the transition operator T in Eq. (10), reading 
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(p|T 



47r 



^(2£+l)P,(p-p') 



E 



1 + mX^fip 



Xa 



2 oo 

E 

n=0 



E{p) - E^e^ 



(20) 
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